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Abstract 

The renormalized expectation value of the stress energy tensor of the confor¬ 
mally invariant massless fields in the Unruh state in the Schwarzschild space- 
time is constructed. It is achieved through solving the conservation equation 
in conformal space and utilizing the regularity conditions in a physical metric. 
The relations of obtained results to existing approximations are analysed. 
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I. INTRODUCTION 


The expectation value of the stress-energy tensor of the conformally invariant massless 
fields in the Unruh state in the Schwarzschild geometry is known to possess some general 
features presented for the first time in the celebrated Christensen and Fulling paper |[J. The 
asymptotic behavior of tangential and radial components of < Tjf > ren and the regularity 
conditions on the future event horizon are quite restrictive, and allow construction of a 
class of approximate tensors. Further, the numerical data, such as the exact value of the 

< Tq > ren on the future event horizon and the value of the luminosity may be used in 
the final determination of the model. It is very fortunate, that thanks to the excellent 
numerical analysis carried out in f2| and || we have detailed informations concerning the 
overall character of the exact < > re n of the scalar field and consequently the validity of 
constructed approximations may by verified. Similar calculations has been carried out in 
the case of the conformal vector field. 

The attempts to construct analytical or semianalytical approximations of < T^ 1 > ren 
in various “vacuum” states, see Refs. [4-19], are motivated, besides self -evident curiosity, 
by the fact that they may be used as a source term of the semiclassical Einstein field 
equations [20-30] or give rise to the analyses of local and averaged energetic conditions and 
quantum inequalities [31-38]. The stress energy tensor is also useful in the thermodynamic 
calculations, see for example [20,25,27,29,39,40]. It should be noted that the back reaction 
calculations are limited to the < T'jf > ren evaluated in the Hartle - Hawking state. 

Recently such semianalytic approximations of the stress tensor in the Unruh state have 
been presented [16-18]. In Refs. [16,17], where properties of the one-loop effective action 
under the conformal transformations have been used, the starting point is the assumption 
that in the ultrastatic companion of the Schwarzschild metric the tangential component of 

< >ren has a polynomial form 

N 

< Tq >ren= Tp(s ) J2 a rU' ? \ (1) 

71—4 
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where x = 2 M/r and T = Tfj /90 (T# = 1/(87 tM)). The numerical factor p(s) depends 
on the spin of the held and is given in the Table 1. Such a choice is in accord with the 
asymptotic analyses of Christensen - Fulling. Further, integrating the conservation equation 
in the conformal space and making use of the regularity conditions in the physical one, a 
family of the stress tensors could be easily obtained. Taking IV = 7 as have been done in 
Refs. [16,17] leaves two free parameters which are to be determined from known value of 
the tangential pressure on the event horizon and the luminosity. It should be noted that it 
is, in a sense, a minimal choice. On the other hand, in Ref. |18| Visser chooses to work in 
the physical space and uses a decomposition of the stress tensor allowed by the covariant 
conservation equation. The tangential component of < TjJ > ren is also taken to be of the 
form (1) with IV = 6, this time however in the Schwarzschild spacetime. The resulting 
model has three free parameters which are determined by performing global unconstrained 
fit to the numerical data. The tensors obtained in both models have a similar structure but 
differ in the numerical coefficients. 

In this note we shall show that, although differently motivated, Visser’s model may be 
obtained following the steps of Refs. [16,17] with N = 8. The simplest N — 5 models 
proposed by Vaz JRJ and by Barrioz and Vaz |TT[ has no free parameters, it should be noted 
however that the latter uses different decomposition of the stress tensor in the conformal 
space. In our constructions scalar, spinor, and vector fields are treated simultaneously. 
Presented method (with n = 1) may be also employed in similar calculations in Hartle - 
Hawking state. 


II. CONFORMAL TRANSFORMATIONS OF THE STRESS - ENERGY TENSOR 


Transformational properties of the stress tensor of a conformally coupled massless held 
under the scaling transformation = e~ 2 ' Ji g IMV are described by or may be obtained from 
the formulas derived by Page |Q, Brown and Ottewill |J Brown, Page, and Ottewill [0, 
and Dowker [T!]]. Under the conformal transformation the one loop effective action of the 
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massless conformally invariant fields transforms according to the rule || 


W R [g^ u ] = W R [e g^] + aA[u\g] + bB[v,g] + cC[u,g], 


( 2 ) 


where 

A[u,g] = 


j d^x^—g) 1 ^ 2 | (^Riem 2 — 2 Rice 2 + ^ [R- + 3 (Do; — «)] (Do; — k) j , 

( 3 ) 


B[uj,g\ = J d A x{—g) 1 / 2 (Riem 2 — ARicc 2 + R 2 ^ ui + AR^u'^u’ 1 ' — 2Rn + 2 k 2 — 4 kUu 


C[u,g\ = J d 4 x(—g) 1 / 2 {[i? + 3 (OR — «)](□/? — «)}, 


( 4 ) 


( 5 ) 


and k = uj- a uj’ a . We have distinguished quantities evaluated in the conformal space with a 
tilde. The spin-dependent coefficients a, b and c are given in Table 1. 

Functionally differentiating (2) with respect to the metric tensor and restricting to the 
Ricci-flat spaces one has 


< > ren = exp (—4a;) Tjf + a(s)A» + b(s)B» + c{s)C^ 


( 6 ) 


where 

A** 


8 R a ^u. aP - + 2g< JV (2ur a K. a + k 2 + - 8- 8uVk, (7) 


B> w = 8R a ^Lu. af3 + 8R a ^uj.. a uj..3 - 8uj^ a u.J u - 8k'^uj' v) - 8 kuj'V 




+ 4 g^ u v. a puj ,a ^ 3 + n. a LJ’ a + - k 2 ^ ’ 


( 8 ) 


and 


C^ = g^i2 Uk + 3 k 2 + 6 uj. a K' a ) - 12 ku'V - 12 k’^lu’A - . ( 9 ) 


The renormalized effective stress tensor has been defined by 


4 




< V 11 ' > ren = 2 (-g) 


( 10 ) 


-1/2 SWr 
tidin' 

We observe that since (6) is a general formula its meaning is clear: the better approx¬ 
imation of the stress tensor in the conformal space is constructed the better < > ren in 
the physical is obtained. There is nothing less general in working in the conformal space 
as long as the transformation formulae are correct. Moreover, although we heavily used 
scaling properties of the renormalized stress tensor the adopted method is neither the Page 
nor Brown and Ottewill approximation. 

Taking the conformal factor in the form u = 1/2 In |r/ 00 |, one obtains for A and 
C'Z the following formulae 


t x 6 (128 - 240 a; + 113 a: 2 ) 

128 M 4 (1 — x) 2 

(ii) 

a; 6 (32 - 96 a;+ 63 a; 2 ) 

381 .1/ 1 (1 - .r) 2 ’ 

( 12 ) 

B -a: 6 (64 - 120 a;+ 57 a; 2 ) 

6 ~ 384 M 4 (1 — x) 2 

(13) 

t 3 x 6 (48 — 96 x + 47 x 2 ) 

‘ ~ 128 M 4 (1 — x) 2 

(14) 

a; 6 (16-48 a;+ 33 a; 2 ) 
r _ 128 M 4 (1 — x) 2 ’ 

(15) 

b q -a; 6 (32 - 72 a:+ 39 a; 2 ) 

9 ~ 128 M 4 (1 — x) 2 

(16) 

t 3a: 6 (32 - 64a; + 33 x 2 ) 

4 _ 256 M 4 (1 — x) 2 

( 17 ) 


— (8 — 9 a;) 2 a; 6 
256 M 4 (1 — x) 2 ’ 


(18) 
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(19) 


0 _ a: 6 (128 — 264 x + 135 x 2 ) 

9 ~ 256 M 4 (1 — x) 2 

The geometrical terms, i. e. the sum of the last three terms in the right hand side of Eq. 

(6) have, therefore, a simple form, proportional to simple polynomials in x multiplied by the 

second power of g tt . 

It is a well known property of the massless and conformally coupled fields that although 
on the classical level the trace of the stress energy tensor identically vanish, after quantization 
it acquires nonzero value. Making use of the identity 

^-S[e~ 2u g^\ ^=0 = —‘2g<jTT-—S\g f j, l ,\, (20) 

ou og aT 

where S is a functional, one concludes that the trace anomaly is given by a general formula 
< > ren — a(s) (^H + — OR^j + b(s)Q + c(s)d.R, (21) 

where 

7 H = Riem 2 — 2 Rice 2 + l-R 2 (22) 

o 

and 


Q = Riem 2 — 4 Rice 2 + R 2 . (23) 

It is a very fortunate coincidence that for the scalar and spinor field < T > ren in the optical 
metric is zero. However, choosing for the vector field the coefficient c(l) as predicted by 
(— function renormalizaton requires some care. Indeed, in this case the trace of < T'jf > ren 
does not vanish and, of course, further calculations should reflect this fact. 

The conservation equation for the line element of the conformal space 

V, l <T^> ren = 0, (24) 


reduces to 


— <r> 

dx < ’ 


2(1 — 2x) 
x(l — x) 


<T[> 


ren 


o, 


(25) 
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and 


— < T r > — ^ - X < T r > + ^ ^ x < T e > — 0 (26) 

dx r ren x(l-x) r ren x(l-x) 6 ren [ ’ 

Therefore one concludes that the stress tensor in the conformal space naturally splits 
into two parts 

dr 

< t; >„„ = V‘ + (27) 

where TJ 1 is a conserved traceless tensor in the conformal space. Such a decomposition is 
allowed by the covariant conservation equation in the ultrastatic space and its exact form is 
constructed from (21) evaluated in the optical space. Similar term has been introduced into 
the Page approximation by Zannias [6]. It should be noted that in a case of the dimensional 
renormalization the second term in the right hand side of the equation (27) is absent. 


III. THE MODEL 


As has been said the idea of constructing < T ^ > ren in the Schwarzschild spacetime is to 
accept (1) and solve the conservation equation in the ultrastatic space. In order to reduce 
the number of unknown parameters we transform (1) to the physical space and make use of 
the regularity condition 


lim I < Tg > ren | < OO. 

X —>1 


(28) 


Inserting (11 - 19) into (6) with appropriate numerical coefficients and making use of (28) 
yields 


a b c 

CLy — ~TZ — — 7T~ — + ——— — 4(24 — 3(25 — 2(26 


48 pt 32 pt 32 pt 


(29) 


and 


7a 7 b 7c 

a s —-T - T - — 7 T 3ct4 T 205 T 

8 384pf 128 pt 256 1 


(30) 
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where t = M 4 T. After substitution of Eqs. (29) and (30) into (1) and performing integration 
of the equations (24) and (25) in the conformal space one has 


< T r > 

^ ± i ^ren 


= — x 2 (l — x) 


, K 
4M 4 ’ 


(31) 


and the radial component of the stress tensor in the conformal space 

a b 


< T r r > ren = pT 


+ 


384 pt 128 pt 256 pt 


— CL4X + 

t b 


-a 4 - -a 5 ) x + 


— (1 I x — 
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+ 


384p t 128 p t 256 p t 


— d\x A 


3° 4 12 a5 ~ 2° 6 ^ X ~ 


16 9 


29 
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^40 pt 
( 7 


120 pt 
21 b 


80 pt 
21 c 


640p t 640 pt 1280p t 


T + 5° 4 + 20° 5 + 10 fl6 J * 

21 9 6 3 \ 8 

___ a4 _ a5 __ a6 u 


(32) 


where K is an integration constant connected to the luminosity and d is another integration 
constant. Note that the leading behavior of (32) as a; —> 0 is proportional to x 2 as expected. 

Now, on general grounds one expects that at large r the leading terms of — < Tj; > ren 
and < 77 >ren should be equal to < > ren , where r* is the Regge - Wheleer coordinate. 
Moreover, the Christensen - Fulling regularity conditions in the Schwarzschild space, i. e. 
conditions for the regularity of the stress - energy tensor on the future event horizon 


< T, m > 7 


< 00, 


(33) 


< I? >ren + < V >ren | < OO, 


(34) 


and 


(1 x) | T uu ^ren | ^ 


(35) 


as x —> 1, where 


Tuu ren 


(1 - x)(< T r r > ren - 


< t; >r 


l) + 


2 K 


(36) 


and 



^ Tw ren ^ T uu > r 


K 

M 2 r 2 


allows to reduce the number of unknown parameters to three. This conditions together with 
(28) ensure that the stress-energy tensor measured in the local frames on the future event 
horizon will be finite. 

The conditions (28) and (34) are already satisfied while the remaining ones after simple 
calculations give 


/ 14 17 1 \ 11a 6 c 

K = pt[--a 4 - -a, - -a e ) 


, 7 17 1 83a 37 b 37c 

Cl — Cl a -\- - ~\~ Cl a — - — - -. (39) 

30 4 120 5 20 6 1920 pt 640 pt 1280 pt v ; 

Returning to the physical space, after some algebra, one obtains the mean value of the energy 
momentum tensor of the quantized, massless, conformally invariant held in the Unruh state 


, pi \ c 11a 6 7 17 a 6 \ 9 , 

<C T. > ren = - \ - — - — - T -a.4 T -a5 T - X — (l 4 X 

4 en l-x\\40 pt 240 pt 20 pt 30 120 20/ 

(5 4 \ 5 /a 9b 3c 7 3 \ 6 

^a i+ -a s )x + [~ t + sft - ^7 - - 2 a y x 

( 7c 109a 416 21 14 7 \ 7 ) 

+ lioy ~ i2o pi ~ Wt + y" 4 + y as + 5° a ) x / ( 


< T r > 

^ - 1 r ^ren 


pT (f 11a 6 c 7 17 a 6 \ 2 4 

1- x\\240pt 20 pt 40 pt 30 120 20/ 

/1 2 \ 5 / a 6 cl l \ 6 

- ( 304 + 305 )* + (iyi + - 4b “ 4 ° 5 - -2“»j X 

/3c 9a 96 9 6 3 \ 7 ) . 

V 10 pt 40 pt 40pt 5 5 5 5 7 j V 


< Tg > ren — pT < a ^ 4 + ( 2 a 4 + a. 5 ) x 5 + -—- 


a 6 \ 

—— — -—- + 3a 4 + 2 a 5 + a@ ] x 
opt 4pt ) 


9 



( 43 ) 


K 

< T r > = —x 2 _ 

^ l t >ren X 4M 4 ’ 


where K is given by 


K = — (11a + 126 - 6 c) 
60 


vT 

(28a4 + 17a 5 — 6 a 6 ) 


(44) 


Generalizations to greater N are obvious, however, it seems that the more complicated 
formulae are of little use. It should stressed that, by construction, obtained tensors satisfy 
all regularity and consistency requirements. 

IV. DISCUSSION 


In order to compare just obtained < T ^ > reri with those constructed by Visser in Ref. (18 
first we introduce a new set of unknown parameters ki defined as 


04 = &4, 


(45) 


05 — — 2k 4 + A? 5 , 


(46) 


and 


1 fa b c' 

o 6 = ki ~ 2h + h + y t u + i~2 


In terms of k , the stress tensor becomes 

„2 ( u u 3 


<■ T 1 * ■> — 

^ t ^ven 


pTx _ 

1 — (c (24 20 80 pt 

3 . 5, 3 

-—- (a + b) + — 

Apt 4 


(a + b) — k±x 2 + (k 4 ~ —k $^ x 3 + 


(a + b) + jk 5 - ^ k e x 4 + 


7 


27 


-ka — 

5 40 pt 


(a + b ) 




(47) 


(48) 


<T r > = 

^ -l r ^ ren 


Ehf (_1_ (0 + 6 ) - h _ h _ ^ + P 4 _ £ fcs j 

1-x [80pt V ; 24 20 V 3 V 


x 3 + 


h 

2 


x 4 + 


AOpt 


('CL + b) 


X 


(49) 


and 
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< Tg > r en= pTx 4 (k A + k 5 X + k 6 X 2 ). 


( 50 ) 


The net flux is described by (43), where the integration constant K is expressed in terms of 
ki as 


K = — [9a + 9b — 2 (5 k 5 + 6 k 6 ) pt\. 
60 


(51) 


Eqs(48 - 51) are equivalents of the < TjJ 1 > ren constructed recently by Visser [18 


The parameters a 4 , a 5 , a 6 or equivalently fc 4 , k 5 , k 6 are to be determined from the avail¬ 
able numerical data. Such a procedure has proven to be very useful in constructing highly 
accurate analytical approximations to the exact stress tensor. In it simplest form one needs 
the horizon value of, say, < Tjj > ren the luminosity and one additional piece of infor¬ 
mation. When the results of detailed numerical calculations are known one may perform 
unconstrained fit to the totality of the available data. In this approach known value of the 
luminosity is used by (51) as a consistency check rather than an input. This procedure has 
been adopted recently by Visser in the case of the scalar field. 

Unfortunately, detailed calculations, even if executed, are rarely published, rather, the 
overall character of the stress tensor is presented graphically. However, in the static and 
spherically - symmetric geometries it is relatively easy to construct the asymptotic charac¬ 
teristics of < > ren • Consequently, with < Tfj( 1) > ren and the luminosity treated as input 
one may construct reasonable approximation. Indeed, taking into account a more restrictive 
hypothesis N — 7, one concludes that the additional constraint 


3a 4 T 2a 5 T ag 


1 / 21 

-( 7 a + 215 -|-c ) = 0, 

384 pt V 2 


(52) 


results, after some rearrangement in the stress tensor 

„2 r , 119 


<* r r t > — 


X 


4M 4 (1 — x) L 


K + 4 h + 24 K - 


32 


123 27 

a — -b —c | ar 

32 64 


303 


339, 


234 


Ah + 56 K — - a — - b + 

32 32 64 


c x 


297 

357, 

405 \ 4 

bll3 

141, 

189 \ ,1 

( 30 A — - a — 

- b + 

c)x 4 + 

- a + 

- b - 

- c I x 5 

V 32 

32 

64 J 

b 32 

32 

64 J 


, ( 53 ) 
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< T r > 

^ - 1 r ^ren 


X 


4M 4 (1 - x) 

211 

+ [Ah + AOK - —a 


K — [Ah + 24 K - —a - —b + — c | x‘ 

32 32 64 


32 


321, 135 , o 

- b H- c x 6 

32 64 ' 


18 AT 


113 141, 189 \ 4 /21 33, 81 , = 

- a — - b + ——c ) x — —a + —b — —c x 

32 32 64 J V32 32 64 


( 54 ) 


and 


x 


4 r 


^ 1 9 ^ren ~ ^ 


, 119 123 27 / 69 81, 81 

h + 6A — - a — - b — -c — 6A — —a — —b H-c | x 

128 128 256 V 64 64 128 


39 

i a - 

V128 128 


/ 19 


b - 


135 

256* 


c x 


(55) 


where h = M 4 (T 9 (2M)). It is interesting to note that both N = 8 and N = 7 hypotheses 
yield the same structure of the renormalized stress tensor in the Unruh state. Substituting 
the spin-dependent coefficients taken from table 1 for scalar and vector fields one easily 
obtains results presented in Refs. [[TbJ and [jT7[| . In a case of the conformal vector field the 
coefficient c(l) has been taken as predicted by C ~function renormalization. Therefore, one 
can draw a conclusion that from the point of view of the applied method the only difference 
between the results of [16,17] and [18] is the choice of N in (1). On the other hand, Barrioz 
and Vaz Jll| take N — 5 ( i. e. there are no free parameters left) and use more complicated 
decomposition of < T> ren in the optical space 


< > ren = 7? + (a« + P5£)< T° 


> r 


(56) 


where a and (3 are coefficients subjected to the obvious condition a + A/3 = 1. 


V. CONCLUDING REMARKS 

In this work our goal was to construct < TJf > ren in the Unruh state and to investi¬ 
gate how the choice of N in Eq. (1) affects the resulting stress-energy tensor. Although 
our analyses have been limited to N < 8 it seems that a three parameter family of the 
stress tensor is of sufficient generality. Since apparently the ambitious plan to construct the 
approximate stress tensor in the Unruh state using the polynomial ansatz and appropriate 
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regularity conditions as the only available a priori informations has failed it seems that 
the presented method (or the methods closely related) are the only one which would give 
analytical formulae able to reconstruct the exact < Tjf > re n to a high accuracy. Moreover, 
the price one should pay for such quality of the approximation is rather small: just two or 
three pieces of numeric data, As far as we know there are neither numerical calculations 
nor asymptotic analyses concerning the vacuum polarization effect of the conformally cou¬ 
pled massless spinor held in the Schwarzschild spacetime and consequently the stress-energy 
tensor cannot be determined completely. We expect however, that the general formulae 
supplemented by additional pieces of numerical data would give a good approximation of 
the exact stress tensor in this case also. 

Finally, we remark that a similar method, with different asymptotics may be used in con¬ 
struction of < Tjj* > ren in the Hartle - Hawking state, specifically, the results of Refs. |L4| may 


be rederived [421. We intend to return to this group of problems in a separate publication. 
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TABLES 


TABLE I. The coefficients a(s), b(s), c(s), and p(s) for fields of helicity s. Fos s = 1 the coeffi¬ 
cient c is nonzero as predicted by point separation and function renormalization. Dimensional 
regularization gives c(l) = 0. 


s 

0 

1/2 

1 

57607r 2 a(s) 

3 

9 

36 

5760tt 2 &0) 

-1 

-11/2 

-62 

57607t 2 c(s) 

0 

0 

-60 

p(s) 

1 

7/4 

2 


14 







REFERENCES 


[1] S. M. Christensen and A. S. Fulling, Phys. Rev. D15, 2088 (1977). 

[2] T. Elser Phys. Lett. A.94, 205 (1983). 

[3] B. P. Jensen, J. G. Me Laughlin, and A. C. Ottewill, Phys. Rev. D43, 4162 (1991). 

[4] D. N. Page, Phys. Rev. D25, 1499 (1982). 

[5] M. R. Brown, and A. C. Ottewill, Phys. Rev. D31, 2514 (1985). 

[6] T. Zannias, Phys. Rev. 30, 1161 (1984). 

[7] M. R. Brown, A. C. Ottewill, and D. N. Page, Phys. Rev. D33, 2840 (1986). 

[8] V. P. Frolov and A. I. Zcl’nikov, Phys. Rev. D35, 3031 (1987). 

[9] I. D. Novikov and V. P. Frolov, Physics of Black Holes , (Nauka, Moscow, 1986) (in 

Russian). 

[10] C. Vaz, Phys. Rev. D39, 1776 (1989). 

[11] F. A. Barrioz and C. Vaz, Phys. Rev. D40, 1340 (1989). 

[12] P. R. Anderson, W. A. Hiscock, and D. A. Samuel, Phys. Rev. Lett. 70, 1739 (1993). 

[13] P. R. Anderson, W. A. Hiscock, and D. A. Samuel, Phys. Rev. D51, 4337 (1995). 

[14] J. Matyjasek, Phys. Rev. D53, 794 (1996); IF UMCS Report 1995, 33. 

[15] J. Matyjasek, Acta. Phys. Polon B27, 2005 (1996). 

[16] J. Matyjasek, Phys. Rev. D55, 809 (1997). 

[17] J. Matyjasek, Class. Quantum Grav. 14, L15 (1997). 

[18] M. Visser, Phys. Rev. D56, 936 (1997). 

[19] J. Matyjasek, Phys. Rev. D57, 7615 (1998). 


15 



[20] J. W. York, Phys. Rev. D31, 775 (1985). 

[21] C. O. Lousto and N. Sanchez, Phys. Lett. B212 411, (1988). 

[22] W. A. Hiscock, Phys. Rev. D37, 2142 (1988). 

[23] S. Massar, Phys. Rev. D52, 5857 (1995). 

[24] D. Hochberg and T. W. Kephart, Phys. Rev. D47, 1465 (1993). 

[25] D. Hochberg, T. W. Kephart and J. W. York, Phys. Rev. D48, 479 (1993). 

[26] D. Hochberg, T. W. Kephart and J. W. York, Phys. Rev. D49, 5257 (1994). 

[27] P. R. Anderson, W. A. Hiscock, J. Whitesell, and J. W. York, Phys. Rev. D50. 6427 
(1994). 

[28] W. A. Hiscock, S. L. Larson, and P. R. Anderson, Phys. Rev. D56, 3571 (1997). 

[29] D. Hochberg and S. V. Snshkov, Phys. Rev. D53, 7094 (1996). 

[30] J. Matyjasek, Acta Phys. Polon. B29, 529 (1998). 

[31] L. H. Ford and T. A. Roman, Phys. Rev. D48, 776 (1993). 

[32] L. H. Ford and T. A. Roman, Phys. Rev. D51, 4277 (1995). 

[33] L. H. Ford and T. A. Roman, Phys. Rev. D53, 1988 (1996). 

[34] M. Visser, Phys. Lett. B349, 443 (1995). 

[35] M. Visser, Lorentzian wormholes - from Einstein to Hawking, American Institute of 
Physics, New York 1995. 

[36] M. Visser, Phys. Rev. D54, 5103 (1996). 

[37] M. Visser, Phys. Rev. D54, 5116 (1996). 

[38] M. Visser, Phys. Rev. D56, 5123 (1996). 


16 



[39] O. B. Zaslavskii, Phys. Lett. A181, 105 (1993). 

[40] O. B. Zaslavskii. Class. Quantum Grav. 13, L23 (1996). 

[41] J. S. Dowker, Phys. Rev. D33, 3150 (1986). 

[42] J. Matyjasek, unpublished. 


17 



